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ABSTRACT 

We revisit the application of Modified Newtonian Dynamics (MOND) to galaxy clus- 
ters. We confront the high quality X-ray data for eight clusters of galaxies observed by 
the XMM-Newton satellite with the predictions of MOND. We obtain a ratio of the 
MOND dynamical mass to the baryonic mass of M m /Mb = 4.94 ± 0.50 in the outer 
parts (i.e r ~ 0.5 R„ir), in the concordance cosmological model where the predicted 
asymptotic ratio, if any baryons are present, is 7.7T^i (at r ^ 0.3 R V i r ). We confirm 
that the MOND paradigm lowers the discrepancy between the binding mass and the 
baryonic mass in clusters by a factor of ~ 1.6 at about half the virial radius. However, 
at thi s radius about 80% of the mass is still missing, and as pointed out bv ISandera 
(2003j), this necessitates a component of dark baryons or neutrinos in the cluster core. 
Concerning the neutrino hypothesis, application of the new data requires a minimum 
neutrino mass of m v > 1.74 ± 0.34 eV to fill this gap. The corresponding 2a lower 
limit of m„ > 1.06 eV is marginally inconsistent with the current constraints from 
the cluster number counts, and from the CMB and large scale structure data. MOND 
must invoke neutrinos to represent the main component that account for the missing 
mass problem in clusters. 
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^ ■ 1 INTRODUCTION 



The missing mass problem in clusters of galaxies arises from 
the comparison of the observed baryonic mass with the ob- 
served dynamic mass. The baryonic mass is mainly due to 
the hot intracluster gas that is well observed in X-rays via 
its free-free emission. The current status of the observed gas 
fraction in clus ters gives a fairly well constrained value of 
about 12% fsee lGreeo et alJl200 ill Allen et alJl200l for in- 
stance). Taking into account the stellar mass, this makes 
the discrepancy between the observed dynamic mass and 
the observed baryonic mass larger than a factor of 7. 

The dark matter (DM hereafter) hypothesis appears to 
provide a seductive explanation of this problem. A new com- 
ponent of non-baryonic matter, insensitive to all interac- 
tions but gravitation, is introduced to fill the gap between 
the baryonic matter and the binding mass. While cosmo- 
logical evidence is accumulating in favour of this scenario 
(see for instance iFreedman et alll200lt ISperget et al.ll2003t 
iTeemark et al J 12004ft ." it is disconcerting that the nature 
of the non-baryonic dark matter is completely unknown. 
Of course there are many candidates of var ying degrees of 
detectability and plausibility (eg, review by iBertone et alJ 
2001) 



As an alternative to dark matter scenarios, iMilgrorr] 
(1983c) proposed a modification of the Newtonian dynam- 
ics effective at galactic and extra-galactic scales. This mod- 
ified Newtonian dynamics (MOND hereafter) has been no- 
tably successful in explaining the discrepancy be tween ro- 
tation and luminosity curves in spiral galaxies jMilgroml 
1983allbh. and claims other phenomenological successes (see 



Sanders fc McOaughl20o'2l for a review). Given that there is 
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now a relativistic, Lorentz invariant generalisation of MOND 
that can incorporate b oth gravitational lensing and cosmol- 
ogy l)BekensteirJl200'3l . it is timely to reexamine one of the 
few admitted failures of MOND. The discrepancy between 
the baryonic mass and the dynamical mass in clusters of 
galaxies is perhaps foremost among the issues that MOND 
has yet to convincingly address. 

The first confrontation of X- ray observations of clusters 
with MOND llOerbal e t alj 1992) e mphasised, despite some 
minor controversy iMilgromlll993t iGerbal et al1ll993ft . the 
difficulties faced by MO ND in pa s sing the clu ster test. The 
problem was revisited bv lSandersl l|l994Lll999ft and ended in 
a remaining discrepancy of a factor of 2-3 bet ween the bary- 
onic a nd the MOND masses. More recently, lAguirre et alJ 
( 2001) discussed observational evidence for three clusters for 
which the observed discrepancy is about 1-5 within 1 Mpc 
and is boosted to a factor ~ 10 within the central 200 kpc, 
further weakening the r eliability of the MOND paradigm. 
However ISandera <l2003ft responded with an update of his 
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earlier work, introducing an added ad hoc dark component 
at the cluster centre to reduce the discrepancy to only a 
factor of 1-3 overall in the cluster. 

Some other tests have also been carried out using grav- 
itational lensing data. They have also pointe d out the diffi- 
culti es faced by M OND at the cluster scale jGavazzill200l 
IClowe et al.ll2004) . 

In this paper, we test the MOND prediction in terms 
of dynamical mass with respect to the observed baryonic 
mass, by basing our work on high quality X-ray data recently 
obtained from observations with the XMM-Newton satellite 
for ten nearby rela xed clusters iPointecouteau et al.ll2005t 
lArnaud et alJl2005Tl . Section 2 presents the context and the 
formalism within which the mass of clusters is derived in the 
MOND case. In section 3 we present the data and the work- 
ing framework we adopt to compute the observed MOND 
masses. The comparison of the MOND mass with the New- 
tonian dynamic mass and the baryonic mass is developed in 
Sec. 4. The results are discussed with respect to previous 
studies in Sec. 5. Our new result is that, in addition to the 
known problem in the cluster cores, there is a significant 
discrepancy in the outer cluster, where an additional ad hoc 
dark component is needed to rescue MOND. We investigate 
the neutrino hypothesis for such a possibility. 

Unless mentioned otherwise, we choose to work in a 
concordance model, using as cosmological parameters: fl — 
1, Q m = 0.3, Ha = 0.7 and Ho = TO km/Mpc/s (referred as 
LCDM). 



2 MOND DYNAMICAL MASSES 

For a spherical system in hydrostatic equilibrium, the den- 
sity and the temperature distributions are connected via the 
equation of hydrostatic equilibrium, so that the dynamical 
mass of a spherical system can be expressed as : 



M d (r) = 



kT r 
G/im p 
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where G is the gravitational constant, m p is the proton mass 
and /J, = 0.609. Mj,(r) is the dynamical mass enclosed within 
a sphere of radius r. Here, n g (r) and T(r) are the density 
and temperature radial distributions. 

In the framework of modified Newtonian dynamics, the 
gravitational acceleration a is linked to the Newtonian ac- 
celeration g as follows: 



af(a/a ) 



(2) 



where ao is assumed to be a universal constant of acceler- 
ation in MOND. Its value as determined from the rotation 
curves of galaxies is about 10 -8 cm/s 2 . In the follow ing, we 
adopt the value of 0.8 X 10~ 8 used bvlSandersI dl999l) (here- 
after S99 - referring to lBegeman et, aL^^9lF ~The transition 
between the Newtonian and the MOND regimes is expressed 
through the functional f(x), x being a/ao, that is also de- 
rived from the application of MOND to the rotation and 
luminosity curve of galaxies: 



(3) 



gravitational acceleration a, quantity that decreases with 
increasing radius: 



a(r) 
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Eq. 8 from S99 relates the dynamical MOND mass, M m , 
to the dynamical Newtonian mass, Md'- 

M d (r) 



M m [r) 



(•») 



Further details o f the MOND formalism can be fo und, for 
instan ce, in S99, ISanders fc McGauehl (120021) and ISandersI 
J2003h . 

According to Eq.^and Eq.^] we can express the ratio 
of the Newtonian dynamical mass to the MOND mass as: 



M d = 
M m (T> 



1+1 ^ — 

G Md 



(6) 



We now quote Eq. 9 from S99, that gives the MOND 



In the following we refer to the Newtonian dynamical mass, 
Md(r), as the dynamical mass, and to the dynamical mass 
in modified Newtonian dynamics, M m , as the MOND mass. 



3 WORKING FRAMEWORK 
3.1 The X-ray data 

Recently, IPointecouteau et all J2005I . hereafter PAP05) 
have published the observed mass profiles for a set of ten 
nearby [z ^ 0.15) and relaxed galaxy clusters observed with 
the XMM-Newton satellite. These clusters, chosen for their 
regular X-ray morphologies, have dynamical masses deter- 
mined through the hypothesis of spherical symmetry and 
the use of Eq. Q They cover a large temperature range of 
[2 — 9] keV, corresponding to a dynamical mass range of 
[1.2 — 12.0] x 1O 14 M0 enclosed within a radius -R200 (i-e a 
radius encompassing 200 times the critical density of the 
Universe at the cluster redshift), and a range of radii be- 
tween 0.01 to 0.7 in units of i?/i?200- 

From the initial sample of ten clusters, we kept only 
eight for the present work: A1991, A2717, A2597, A1068, 
A478, A1413, PKS0745, A2204. We excluded A1983 and 
MKW9. First they are the only two clusters of the sample 
that are not observed up to 0.5 i?200- Indeed, we need to keep 
to the observed radial range of each cluster to avoid any ex- 
trapolation, and to derive reliable observational constraints. 
Meanwhile MWK9 has also the most disturb ed morphol- 
ogy of the 10 clusters jPratt fc Arnaudl 12005^ . Both also 
exhibit unexpectedly low gas fractions that may turn them 
into outliers in terms of the average gas fraction in clusters 
\Greeo et all200l1 : lAllen et all2003t) . As our sample is quite 
small, excluding those two clusters may avoid any poissonian 
bias on the computation of the intrinsic dispersion for the 
ratio of the MOND mass to the baryonic mass. However, 
at the scale of our sample, this also may slightly bias the 
results in favour of MOND predictions, for which low gas 
fraction systems will be difficult to explain (as it increases 
the discrepancy between the observed baryonic mass and 
the dynamical mass). 

We used the observed density and temperature profiles 
derived by PAP05 to compute the observed ratios of the 
dynamical mass to the MOND mass and of the MOND mass 
to the baryonic mass. 
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Figure 1. Observed scaled profiles for the ratio of the dynamical 
mass to the MOND mass. The average profile over the sample 
is overplotted (solid black line) together with its lcr dispersion 
(dashed lines). The triangles mark the estimate from the self- 
similar hypothesis for each cluster (see Sect. 

3.2 Scaled MOND profiles 

In hierarchical structure formation, the virialized part of 
clusters at a given redshift is encompassed within a sphere 
with a radius corresponding to a fixed density contrast 
(~ 200) with respect to the critical density of the Universe. 
-R20 is thus considered to be the virial radius of the cluster 
(see lArnaud et al]l2004l and references therein). 

Thus to cross-compare the internal structure of clus- 
ters, one has to look at quantities scaled according to ifeoo- 
As an extension, to compare the MOND and the dynami- 
cal masses as a function of radius, we choose to scale both 
quantities according to the virial radius (i.e Jfeoo)- For each 
cluster, we used the value of R200 , reported by PAP05 for a 
LCDM cosmol ogy (see their table 2 ), obtained from the fit of 
NFW profiles avarro et al]ll997l) over the observed mass 
profiles. In the f r amew ork of MOND structure formation, 
iKnebe fc Gibson! (120041) have shown that the most massive 
galactic halos formed in their numerical simulations could 
be fitted by a NFW profile. We could thus fit each MOND 
mass profile with a NFW profile to derive the correspond- 
ing characteristic scaling radius -R^oo- However, for a given 
cluster, the values of -R200 and f?|oo are likely to be differ- 
ent. Scaling Md(r) and M m (r) respectively with -R200 and 
R200 will lead us to compare these two profiles at different 
physical scales. For our purpose, working with scaled radii 
is required by the use of a sample of clusters with different 
masses, but is not required to compare the MOND and dy- 
namical masses for a given cluster. A single scaling has to 
be adopted for each cluster.We thus adopt -R200 as a scaling 
radius for both dynamical and MOND profiles. 

3.3 The baryonic mass in clusters 

The observational constraints on the gas fraction in clus- 
ters strongly favo ur an asymptotic value of about 10- 12% 
beyond 0.3 R200 llGreeo et al.ll200ll : lAUen et all 120031) . in 
a LCDM cosmology. For the sample used here, the aver- 
age gas fraction at a density contrast of 8 = 1000 (i.e at 



r = 0.47±0.02 R200 - see APP05) is < f gaB >= 0.11±0.02. 
We adopted this value in the following, as the average 
asymptotic value for the gas fraction. 

We now can further investigate the issue of the stel- 
lar mass. For the stellar mas s, S99 used a corre lation de- 
rived from the early work of iDavid et alJ l|l990F ) between 
the luminosity and the gas masses in clusters leading to: 
/* = M^/Mgas ~ 0.7(fcT/keV)~ 1 /i 3/2 . We used this esti- 
mate when working in CDM cosmology, making use of the 
spectroscopic temperature measured between [0.1, 0.5] R200 
(see APP05). Nevertheless, more recent wo rk allows u s to re- 
vise this crude estimate of /*. For example. iLin et, ail (|2003ll 
derived a variation of /* with the total mass, and thus with 
the temperature. From their equation 10 and the values of 
M500 derived by APP05 (see their table 1), we obtain an 
average value of < /* >= 0. 12 ± 0.03 (values ranging from 
0.1 to 0.17). In related work. IVoevodkin fc Vikhlininl feooih 
obtained consistent, but slightly higher, values, and con- 
cluded that the stellar mass represents about 15% of the 
gas mass. We adopted this value of /* = 15% for our sam- 
ple. The baryon fraction will then be /& = 0.13 ± 0.02. In 
order to cope with the baryonic mass in clusters, the MOND 
mass has then to be 7.Tt\'o times lower than the dynamical 
mass once the baryonic mass reaches its asymptotic value 
(i.e r > 0.3 R200). 



4 MATCHING THE BARYON CONTENTS OF 
CLUSTERS 

In an initial approach, taking into account the self-similar 
nature of the cluster population, we can express the ratio 
MdjMm as a function of kT. We make use of the scaling rela- 
tions between the mass and the temperature: M = A(z) Tg" , 
and between the mass and the radius M = iiv/3 S p(z) R\. 
Then Eq. El can be rewritten as follows: 



= Vl + BWTf- (7) 

where B(z) = {3ao / AnS p c (z)) A(z)~ , and Ts is the tem- 
perature at density contrast 5. 

Using the R~T relation published by APP05 at 5 = 
1000 (a density contrast corresponding to observed radii for 
all eight clusters) we estimated from the above equation 
M d /Mm = 1.63 ± 0.29, which leads to M m /M b ~ 4.7. 

We can reverse the computation and estimate the 
value of ao needed to reach an equality between the 
MOND mass M m and baryonic mass Mb. We obtain 
a = (4.75 ± 1.24) x 10~ 8 cm/s 2 at S = 1000 using the eight 
clusters. 

From the mass profiles derived by PAP05 in a LCDM 
cosmology, we computed the observed ratio Md/Mm from 
Eq.|S]for the eight clusters of the sample. The profiles (scaled 
to R200) are shown in Fig. The estimates from the self- 
similar hypothesis are shown for <5 = 1000 as the observed 
triangles, and match the computed ratio very well. 

With the previously given notations and Eq.|S| we can 
express the ratio of the MOND mass to the baryonic mass 
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Table 1. Ratios of the dynamical mass to the MOND mass and 
of the MOND mass to the baryonic mass. 



(8) 



In Tabled we report the average baryon fraction values and 
the average M m /Mb and Md/M m ratios for the measured 
clusters at a given radius. The computation were done at 
the following radii: the physical radius of 750 kpc in order 
to directly compare with S99, and the radii corresponding to 
the density contrasts of S = 1000 and 5 = 15000 (i.e 0.47 ± 
0.02 -R200 and 0.10 ± 0.01 R200 average over the sample). 
Those two radii mark the boundaries of the radial range 
over which the observational constraints are especially well 
tied down (see PAP05, APP05). 

From Fig. it is obvious that the ratio Md/M m de- 
creases, thus the discrepancy between the MOND and the 
baryonic mass increases, towards the centre of clusters as 
the acceleration returns to the Newtonian regime. At ~ 
0.1 R200, the ratio is 1.08 ± 0.06. The discrepancy also in- 
creases for massive (i.e hot) systems. Indeed if we only con- 
sider the 5 hot clusters 1 we obtain at 8 = 1000 a ratio of 
Md/M m = 1.43 ± 0.08, an even smaller value. The corre- 
sponding value of ao to match the baryonic mass has to be 
a = (5.58 ± 0.63) x 10" 8 cm/s 2 . 



5 DISCUSSION 

5.1 MOND as a stand alone solution 

For a direct comparison with S99, it is interesting to per- 
formed the computation at the radius of 750 kpc with fi m = 
1 and Ho = 50 km/s/Mpc. For the five clusters observed be- 
yond 750 kpc (i.e. the 5 hot clusters), the average gas frac- 
tion is then 0.19 ± 0.02 and the derived ratio of the MOND 
mass to the baryonic mass is M m /Mf, = 3.84 ±0.33. For the 
same cosmological setup, S99 claimed that the MOND con- 
text reduces the discrepancy with the baryonic mass down 
to an "acceptable" factor of ~ 2. However, our value of 3.8 
is significantly higher, and does not even agree with a factor 
of 2 within a 3cr limit. 

As we are working with a sample of nearby clusters, 
the most important cosmological parameter is the Hubble 
constant. Indeed, the baryon fraction in clusters scales with 
h~ 3 / 2 , thus decreases with increasing values of Ho- Current 
observational evidence strongly favours a flat Universe with 
a low matter density Q m ~ 0.3, a component due to dark 
energy which can be represented by a positive cosmological 
constant A, and a hig h value for the Hubble constant (Hp ~ 
70 km/s/Mpc) (see iFreedman et alj |200ll : [Spergel et alJ 
120031: iTegmark et all 120041 : ISeliak fc et all 120051 for in- 
stance). In such a universe (i.e. LCDM) Md/M m decreases 
to 1.63 ± 0.31 (1.08 ± 0.06) at ~ 0.5 R200 (~ 0.1 R200), 
so that the discrepancy is increased between the MOND 
dynamic mass and the baryonic mass. The corresponding 
ratio of the MOND mass to the baryonic mass is now 
4.94 ± 0.50 (10.6 ± 3.77). This is more than a factor of 



1 i.e A1413, A478, A2204, PKS0745, A1068 
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(a) Number of clusters used, (b) Ratio for the hot clusters only 
(i.e kT > 3.5). 



two above the value derived by S99. The evidence is con- 
firmed if we only consider the hot systems. Indeed, for clus- 
ters with (kT > 3.5 keV), we derive M d /M m = 1.43 ± 0.08 
at ~ 0.5 R200, which makes the ratio of the MOND mass to 
the baryonic mass ~ 5.10 ± 0.56. Thus in all cases, within 
a 3<r (i.e. 99% confidence) the ratio M m /M b in a MOND 
cosmology will be greater than 3.4, making MOND unable 
to fully overcome the missing mass problem in clusters. 

The case of A1413 is even more eloquent. Indeed, this 
cluster has been observed beyond -R500 (i.e. a physical radius 
of 11 29 kpc corresponding to ~ 0.7 R200 - iPratt fe Arnaudl 
2002), a radius that conservatively is often taken as the outer 
bound of the virialized part in clusters. Moreover, for this 
cluster the XMM-Newton results agrees very well with the 
results derived from the Chandra observatio ns, especially in 
term of the shape of the temperature profile l|Vikhlinin et alJ 
2005). At the radius of -R500, the measured gas fraction 
is f g as = 0.15 ± 0.01, which makes the baryon fraction 
0.17±0.01 (see Sect. 3.3). We derived a corresponding ratio 
Mm/Mi, = 3.6 ± 0.7. Such an observed discrepancy mea- 
sured within a radius closing the virial radius puts a very 
tight constraint on the MOND paradigm. 

In fact, S99 relies on hypotheses that together might 
have biased its results down towards an optimistic value: (i) 
the use of a /3-model to describe the gas distribution down 
to 750 kpc contributes to overestimate the gas mass at large 
radii. Moreover the spatial resolution of the current X-ray 
data has ruled out the /3-model as a fitted representation 
of the observed X-ray surface brightness profiles (especially 
in nearby clusters - see PAP05 for instance); (ii) The use 
of physical radius as a working radius could have induced 
biases as clusters of very different masses are then compared 
at various different density contrasts; (iii) The isothermality 
of the intra-cluster medium was forced by the single over- 
all temperature measurements available for the considered 
sample. However, in contrast we have made use of accurately 
measured temperature profiles. 

If we consider the problem in terms of the value of the 
characteristic MOND acceleration, we derive values of 00 
that are 5 times larger (and still 3.5 times larger within the 
3cr limit) than the value derived from the rotation curves of 
galaxies, and are thus unacceptable values. 



5.2 Adding a dark component 

A last alternative to rescue MOND is to invoke a non- 
luminous c omponen t at th e centre of clusters, as suggested 
recently bv lSanderj (l2003h (S03 hereafter). This author in- 
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traduced this dark component on an ad hoc basis to ex- 
plai n the huge observed discrepancy within the cluster cen- 
tre l|Aguirre et al.l200lf) . Indeed, the baryon fraction is lower 
in the inner parts, so that the ratio to the MOND mass in- 
creases. With respect to our sample, if such a component 
was to exist it would have to account for 81% (±4%) of the 
total mass at ~ 0.5 R,2oo- In other terms, this means that 
MOND just reduces the missing mass problem in clusters 
by about 20% (at half the virial radius) compare to a stan- 
dard DM scenario, but does not solve it. Therefore, some 
exotic dark component has to be added to fill the remaining 
gap of 80%. This component will become the main explana- 
tion for the missing mass, and will thus draw a cosmological 
setup that will be closer to a mixed DM scenario than to a 
dominating MONDian scenario. 

It is highly unlikely that 80% of the missing mass in 
cluster can be due to hidden baryons. Indeed, if atomic or 
molecular hydrogen exist within the intra-cluster medium, 
it will have to face its thermal conditions. Such a compo- 
nent could only account for a small fraction of the thermal 
baryons. The relativistic populations are also known to be a 
minor component in terms of mass, as otherwise they will be 
expected to produce a strong hard X-ray signal, and stronger 
radio emissions than those currently observed. 

Nevertheless, we can follow the suggestion by S03 to 
consider, as potential candidates for this dark component, 
massive neutrinos that aggregate at cluster scale. In the 
MOND case, according to S03, the needed neutrino mass 
per neutrino flavour needed for MOND to cope with the 
observed baryonic mass is about 2 eV. Further assuming a 
constant density sphere for his dark component and tak- 
ing into account the phase space density limit for neu- 
trinos, S03 derived an upper limit for the neutrino den- 
sity after their collapse and accretion within structure of: 
p v «S (4.8 x 10- 24 )(m„/2eV) 4 (T kcV ) 3/2 kg m" 3 . We make 
use of this limit on the neutrino density, and we use the 
spectroscopic temperatures measured for each cluster of our 
sample between 0.1 and 0.5 -R200 (see PAP05). From our 
eight clusters, it is possible to compute the needed neu- 
trino mass to equate the missing mass at a given radius 
(i.e. M m (r) — Mt(r)) with the contribution of the massive 
neutrino to the cluster total mass. S03 hypothesis of the neu- 
trino accretion mainly concerned the central parts of clus- 
ters. In our study, to explain the ~ 80% of missing mass in 
MOND, we add to extend the radius of the neutrino sphere 
down to -Riooo- The minimum neutrino mass then required 
is m„ > 1.74 ±0.34 eV . This is a strongly constraining value 
for the neutrino mass, which makes the lower bound for the 
neutrino mass becomes ~ 1.06 eV, within a 2a limit (i.e. 
95% confidence). 

To date particle physics experiments (single and dou- 
ble beta-decay, neutrino oscillations measurements) lead to 
a wide range of upper limi t for the neutrin o mass going 
from ~ 0.8 eV to ~ 3 eV iFogli et aljl2004t) . Those con- 
straints directly apply to MOND, and they are consistent 
with our lower limit taking into account their large varia- 
tion. If to date the constraints derived from astrophysical 
data seems to be tighter, they may be considered as irrele- 
vant in a MOND framework (i.e. Sl„/i 2 = m Vi /94eV) as 
they are derived within the framework of a standard cos- 
mological DM model for structure formations. Nevertheless, 
within the cosmological MOND+neutrino setup we end up 



with, MOND plays a minor role, where the neutrinos are the 
major component to explain the missing mass in the Uni- 
verse. We can thus reconsider the astrophysical constraints 
in this context. As the number of clusters is linked to the 
matter content of the Universe, thus to Q m and £2„, the clus- 
ter number cou nts can be used to give an upper limit on the 
neutr ino mass jKahniashvili et alJ 120051 : lElgarOv fc Lahavl 



2005). With the following cosmological setup: Q m = 0.3, 
h — 0.7, n — 1 and 0.7 < as < 1.1 (to accou nt for the varia- 
tion in erg determinations - see for instance iTegmark et alJ 
( 20041: llBeliak et alJ <l2005l) L current cluster number counts 
lead t o an 2a upper lim it of m v < 0.8 eV llFukugita et alJ 
l200(j lAllen et al.tl2003l) . A value marginally incompatible 
with the lower bound was derived here. This inconsistency 
is even stronger with the 2a upper limit of ~ 0.3 eV 
derived from th e combination of the WMAP, 2dFGRS 
and SPSS data JSpergel et al-lEooll; ITegmark et alJliool 
ISeliak fc etalll200Et lElgardv fc Lahavil2005h . However, one 
has to keep in mind that there is quite a large dis persion be- 
tween those upper limits llElgarOv fc Lahavl2005l) . that may 
just leave enough room for the MOND+neutrino hypothesis 
to stand as a reliable paradigm. 



6 CONCLUSION 

We have revisited the case of the galaxy cluster scale in the 
framework of the MOND paradigm. We based our study on a 
set of consistent and recent high quality X-ray data obtained 
with the last generation of X-ray satellites (XMM-Newton ). 
For eight nearby and relaxed clusters, we confirmed that 
MOND alone is not able to explain the missing mass problem 
in clusters of galaxies. Almost 80% of the mass is still missing 
at half the virial radius of clusters in a MOND Universe. 
Thus MOND is not the main solution to the missing mass 
problem at cluster scales. 

Indeed, it undeniably requires extra exotic DM to sur- 
vive. In this context we have investigated the hypothesis of 
massive neutrinos distributed as a sphere of constant den- 
sity up to half the virial radius of galaxy clusters. Under 
this hypothesis, we derived a very tight observational lower 
bound for the neutrino mass, m v > 1.06 eV (95% confi- 
dence) , which is marginally inconsistent with the constraints 
from the cluster number counts and from the CMB and LSS 
constraints. The consistency with the current direct experi- 
ments to measure the neutrino mass is spoiled by their as- 
sociated large uncertainties. The perspectives of dedicated 
CMB and Large-scale structure experiments in coming years 
will bring definitive constraints on m„ in the astrophysical 
context, and should definitively settle the case of the MOND 
paradigm. 

Note that recently ISkordis et alJ i2005t) have been 
studying the formation of structure in the relativis- 
tic MOND frame work (i.e. the Bekenstein theory - see 
iBekensteinl (|2004)). To reproduce the observed angular 
power spectrum of the CMB, those authors appeal to mas- 
sive neutrinos with ra„ ~ 2 eV. So unless the main frac- 
tion of the baryonic content of clusters remains hidden from 
the current observations, in any case a large neutrino mass 
(or another exotic massive candidate) is needed to promote 
MOND as a reliable paradigm. As such, massive neutrinos 
would then become a major component of the Universe mat- 
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ter content. This will turn the cosmological framework more 
into a mixed DM cosmology than into a MONDian cosmol- 
ogy- 
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